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Abstract 

There are two types i = 1, 2 of particles on the line R, with Ni particles of type i . 
Each particle of type i moves with constant velocity Vi. Moreover, any particle of type 
Oh i = 1, 2 jumps to any particle of type j — 1,2 with rates 2V7 ay. We discuss in details the 

^ (J-i initial desynchronization of this particle system, namely, we are interested in behaviour of 

the process when the total number of particles Ni + N2 tends to infinity, N1/N2 — > const 
and the time t > is fixed. 
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1 The Model 



The simplest formulation of the model, we consider here, is in terms of the particle system. 
On the real line there are N\ particles of type 1 and N2 particles of type 2, N = N% + N2. 
Each particle of type i = 1, 2 performs two independent movements. First of all, it moves 
with constant speed Vi in the positive direction. We assume further that Vi are constant and 
different, thus we can assume without loss of generality that < v\ < t>2- The degenerate 
case v\ = V2 is different and will be considered separately. 

Secondly, at any time interval [t, t + dt] each particle of type i independently of the others 
with probability aijdt decides to make a jump to some particle of type j and chooses the 
coordinate of the J-type particle, where to jump, among the particles of type j, with probability 
■jy-. Here a^- are given nonnegative parameters for i,j = 1,2. Further on, unless otherwise 
stated, we assume that an = 022 = 0, ai2, 021 > 0. 

After such instantaneous jump the particle of type i continues the movement with the same 
velocity Vi. This defines continuous time Markov chain i = 1,2; k = l,...,Aj. 

where x^\t) is the coordinate of fc-th particle of type i at time t. We assume that the initial 

coordinates x, * (0) of the particles at time are given. We are interested in the long time 
evolution of this system on various scales with iV — ¥ 00, t = t{N) — > 00. 

In different terms, this can be interpreted as the time synchronization problem. In general, 
time synchronization problem can be presented as follows. There are N systems (processors, 
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units, persons etc.) There is an absolute (physical) time t, but each processor j fulfills a homo- 
geneous job in its own proper time tj = Vjt,Vj > 0. Proper time is measured by the amount 
Vj of the job, accomplished by the processor for the unit of the physical time, if it is disjoint 
from other processors. However, there is a communication between each pair of processors, 
which should lead to drastic change of their proper times. In our case the coordinates x¥\t) 
can be interpreted as the modified proper times of the particles-processors, the nonmodified 
proper time being xj^(0) + Vit. 

There can be many variants of exact formulation of such problem, see [H [JJ [TU| [2J [T5]. We 
will call the model considered here the basic model, because there are no restrictions on the 
jump process. Many other problems include such restrictions, for example, only jumps to the 
left are allowed. Due to absence of restrictions, this problem, as we will see below, is a "linear 
problem" in the sense that after scalings it leads to linear equations. In despite of this it has 
nontrivial behaviour, one sees different picture on different time scales. 

There are, however, other interesting interpretations of this model, related to psychology, 
biology and physics; For example, in social psychology perception of time and life tempo 
strongly depends on the social contacts and intercourse. We will not enter the details here. 



2 Main results 

We show that the process consists of three consecutive stages: initial desynchronization up to 
the critical scale, critical slow down of desynchronization and final stabilization. 



Final stabilization The first theorem shows that for Ni fixed and t — > oo there is a syn- 
chronization: all particles asymptotically, as t — > oo, move with the same constant velocity v, 
that is like vt. However it does not say how fluctuations depend on iVj. 
Put 

m{t) = min xi (t) 

i,k 

Theorem 1 For any fixed Ni,N2 there exists v = v(Ni,N2) > such that for any i = 1,2 
and any k = 1, Ni a.s. 

urn — = v 

t— >oo t 

Moreover, the distribution of the vector ^x^\t) — m(t), i = 1, 2; k = 1, . . . , iVj \ tends to a 
stationary distribution. 

The velocity v will be written down explicitely in terms of this distribution, it depends of 
course on aij and Vi. Note that both the velocity and the distribution do not depend on the 
initial coordinates. 



Initial desynchronization Now we consider the case when N — > oo but t is fixed. More 
exactly, we consider a sequence of pairs (N\,N2) such that N%,N2 — > oo so that -j4- — > Cj, 
where c\ + C2 = l,Cj > 0. It is convenient here to consider positive measures or generalized 
functions 

m^)(i,x) = -^^<5(x-4 i) (i)), xGR + 

k 
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defined by the coordinates of iV; particles of type i at time t. We assume that at time t = 
for any bounded C 1 -functions <f>i{x) on R the sequence < m^ Nl \o, > converges to some 
number. 

Theorem 2 Then for any t there are weak deterministic limits 

lim — m\ Ni \t,x) = rrii(t,x) 

iV— >oo iv 

where rrii(t,x) satisfy the following equations 

dmx dmi 

~df 1 ~dx~ = ai2 ^ m2 ~ ( > 
dm 2 dm 2 

+ v 2 -^r = a 2 i(m 1 - m 2 ) (2) 
Now we want to study the asymptotic behaviour of mi(x,t) for t — >■ oo. Denote 
cii(t) = J xirii(x,t)dx, di(t) = J (x — ai(t)) 2 mi(t, x) dx 
Theorem 3 There exist constants v, d > such that as t — > oo 

CLi{t) =Vt + Oj + 0(1) 

di(t) = dt + djo + o(l) 
/or some constants Ojo,djo- Moreover, 

A l (x,t) = ^4^ 1 

tends to --=exp(— pointwise as t — > oo. 

Critical point and uniform estimates Here we assume that iVi = [ciJV], N 2 = [c 2 N] for 
some Cj > 0, ci + c 2 = 1. Introduce the empirical means (mass centres) for types 1 and 2 

= ^ E 4° (*)> 
* fe=i 

the empirical variances 

N, 



k=l 

and their means 



Sf(t) = ^E{4\t)-x^(t)j 



/x i (t) = ErW(t), Ii 2 (t) = Aii(t)-A*2(t). Ri(t) = ES?(t) 

The following asymptotic results hold for any sequence of pairs (iV,i) with iV — >■ oo and 
f = t(N) -»• oo. 
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Theorem 4 We have the following asymptotical results as t — )■ oo; 

Ul-U 2 £lj(t) 012^2 + 021^1 

hz\t) -)• ■ , — ► ■ 

«12 + «21 t a\2 + OL2X 

Assume now that iVj = CjiV, where Cj > 0, c\ + C2 = 1. 

Theorem 5 There are the following three regions of asymptotic behaviour, uniform in t(N) 
for sufficiently large N: 

• if t -^p- -> 0, i/ien Ri(t(N)) ~ hx 2 t(N), 

• ift = t(N) = sN for some s > 0, &en Ri(t(N)) ~ /i (1 - e^ 2S )iV, 

• i/ -> 00, i/ien Ri(t(N)) ~ WV, 

where the constant k 2 > can 6e explicitely calculated and 

_ 2a 12 a2i {v\ - V2) 2 
^2(^12 + «2l) 3 

The proofs of Theorems [4] and [5] are given in |13[ [14] . 

3 Limit £ — > oo 

In this section we will prove Theorem [TJ 

Two particles. It is useful to consider first the case when iVi = N 2 = 1. Thus consider the 
process (xW(i),x( 2 )(i)). We will prove that there exist deterministic limits 

, x®(t) 
hm = v 

t— >oo t 

for i = l,2 and some v > 0, moreover the distribution of the random variable p(t) = x^ 2 \t) — 
x^(t) tends to some distribution on R+. 

We can assume that v\ = 0, t>2 > 0. The Markov chain p(t) = x^ 2 \t) — x^'(t) on R + 
satisfies the Doeblin condition, that is from any x G R+ there is a jump rate to 0, bounded 
away from zero, here it equals Q12 + 021- It follows that p(t) is ergodic. Then as t — > 00 there 
exists the limiting (invariant) distribution F(x) for p(t). Let 

h <t 2 < ■■■ 

time moments when x^\t) = x^ 2 \t). It is clear that t% — tk-\ are independent random vari- 
ables, exponentially distributed with parameter a 12 + a 2 i- It follows that F(x) is exponential 
with the density 

, x x / \ \ , «12 + a 2 i 

p(x) = Aexpf — Ax), A = 

v 2 - V\ 
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Thus, if the limits lim ^ili exist, then they are equal. Let us prove that they exist and 

t— >oo 1 

v = v\ + ai2 J xp{x)dx (3) 

In fact, the particle 1 moves with constant speed v\ and performs on the time interval [0, T] 
independent exponentially distributed jumps in the positive direction. As T — > oo, the number 
of these jumps asymptotically equals ct\ 2 T, and the mean jump asymptotically is j xp(x)dx. 
Similarly one can get 

v = V2 — a 2 \ J xp(x)dx (4) 

From this and ([H]) we have 

a 2 \Vi + ai 2 v 2 
v = 

«21 + «12 



General case. Let us prove first the second statement of the theorem. We can put V\ = 
and change the coordinate system putting m(t) = 0. Consider a configuration of particles 
at time t. Denote the particle, which has coordinate m(t) = at time t, as particle 0. Let 
p(t + 2) be the probability that at time t + 2 each particle will be inside the interval [0, 2^2]. 
This probability can be (very roughly) estimated from below as 

p(t + 2) > mm(p 01 p 2 pi,po2P3Pi) 

To prove this consider first the case when the particle has type 1. Under this condition 
p(t + 2) can be estimated from below as poip 2 pi, where poi is the probability that particle 
does not do any jumps in the time interval (t,t + 2), p 2 is the probability that each particle 
of type 2 jumps at least once to the particle in the time interval (t,t + 1) and does not 
do any more jumps in the time interval (t,t + 2), p\ is the probability that each particle of 
type 1 jumps to some particle of type 2 in the time interval (t + l,t + 2). Similarly, under 
the condition that the particle has type 2, p{t + 2) can be estimated from below as P02P3P4, 
where pq 2 is the probability that the particle does not do any jumps in the time interval 
(t, t + 2), P3 is the probability that each particle of type 1 jumps at least once to the particle 
in the time interval (t, t + 1) and does not do any more jumps in the time interval (t, t + 2), 
P4 is the probability that each particle of type 2 jumps to some particle of type 1 in the time 
interval (t + l,t + 2). 

This means that the Markov chain C = — m(t),i = 1, 2;k = 1, iVj| satisfies the 

Doeblin condition. Then it is ergodic and has some stationary distribution. We will write 
now formula for v, assuming however that an = 0. For this we need some marginals of this 
stationary distribution. 

Let Ai(t) be the event that at time t at the point m(t) there is a particle of type i, and 

qi = lim P{Ai{t)) be the stationary (limiting) probability of A{. Let pt(y) be the stationary 

t— >oo 

conditional (under the condition Ai) probability density of the distance from m to the nearest 
particle. In the time interval [T,T + dt] the particle in m(t) moves with the speed Vi, and 
moreover can make one jump. This gives, for example under the condition A\, constant 
movement v\dt of m, and the jump of m to the nearest point with rate a\ 2 dt. Thus as 
T — > 00 we have 

E(m(T + dt) I m(T) ) - m(T) = q x ( Wl + a 12 j ypi(y)dy) dt + 
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<?2 (y2 + q.21 j yp2(y)dy^j dt + o(l) 
qi (vi + an J ypi(y)dy^j + q 2 (v 2 + a 2 \ j yp2(y)dy^j 



and then 



About Doeblin chains. In the standard theory of Doeblin chains, see [3], it is assumed 
that transition probabilities are absolutely continuous with respect to some positive measure 
H on the state space. 

(i) (i) 
If at time all x) are different, then for any t it is true that all different a.s. 

Thus transition probabilities (for example, for the embedded chain at times 0,1,2,...) are 

absolutely continuous with respect to Lebesgue measure on (R^ 1_1 x R^ 2 ) U (R+ 1 x R^ 2 ^ 1 ). 

If at time some coordinates coincide, then a.s. in finite time t they become all different. 



4 Limit N 



oo 



It is very intuitive to introduce the following continuous model. Let nii(0,x), x E R, i = 
1,2, be positive smooth functions, M, = J mj(0, x)dx = 1. We call them continuous mass 
distributions of type i at time t = 0. The dynamics of the masses is deterministic — during 
time dt from each element dm\ of the mass the part a\2 dt dm\ goes out and distributes 
correspondingly to the mass 7712(2;), namely it becomes the mass distribution with density 
7712(2;) cni2 dtdmi, and vice- versa, interchanging 1 and 2. Moreover each mass element moves 
with velocities v\ and v 2 correspondingly. From this we easily get linear equations (JI])-© for 
mass distribution rrii(t,x) at time t with the initial conditions 

7774(0,2;) = fi(x) 

Now we will prove convergence of N particle model to the continuous model. 



4.1 Convergence: the martingale problem 

Here we prove Theorem (2) 

We consider continuous time Markov process 

tnM*) = • • .,*$(*); 4 2) (*), • • -,4l(t)) 

with the state space R-^ 1 " 1 "^ 2 . Its generator 



(L Nl>N2 f) (x^-x^ 



EO \ -> 



~^ dx\ dx 



(2) 



f(xW; X W) + 



Ni N 2 



i=l j=l 
N 2 Ni 



iV2 ivi r / v 

.7=1 i=l L V J/ 



a-W. x (2) 



(5) 



+ 
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where the following notation is used 

( r W. r (2)\ _ ( (1) _(1)._(2) (2)\ 

f T (l). r (2)\ _ / r (l) _(1) _(2) (1) r (l). r (2) (2)\ 

^x ,x ~~ ^ x l ' • • • ' x j-l> x j ' x j+l> • • • ' x iVi' x l '•••> x /V 2y /> 

( T (1). T (2)\ (J» r (l). T (2) (2) (1) (2) (2)\ 

^x , x j ^ ^x^ , . . . , x^ '1 ) • • • ' i ' j'+l' ' ' ' ' N2 J ' 

is defined on bounded C 1 -functions. 

We will consider the limiting behaviour of this process when t = const, N±,N2 —> 00. It is 
not convenient to deal with the sequence (,N 1 ,N 2 (t) of processes because the dimension of the 
state space changes with N\,N2- 

Denote 



M Nl:N2 {t) = 




where 6(x), x £ R, is the (^-function. One can see that the generalized functions 

1 i=i 2 3=1 

represent empirical "densities" or masses of (type 1 and 2 correspondingly) particles at 
time t. Thus, if <f>(x) = (fa(x) , fa(x)) , where fa £ S(R), then for fixed particle positions 
Xi(t), . . . ,Xpj(t) and x^\t), . . . ,xffl(t) the vector function Mjy lt N 2 (t) is a linear functional 
on the vector test functions </>, that is 




Fix some T > 0. Then {Msx N 2 (J>)i < t < T) can be considered as a Markov process 
taking its values in the space of tempered distributions S"(R) xS'(R). In the sequel we consider 
S"(R) x 5"(R) as a topological space equipped with the strong topology (see Subsection I5.2p . 
Without loss of generality one can assume that the trajectories of the process Mjv 1 ,Ar 2 (t) are 
right continuous functions with left limits. So it is natural to consider the Skorohod space 
n T = D([0,T],S'(R) x S'(R)) of functions on [0,T] with values in 5'(R) x 5'(R) as a 
coordinate space of the process Mjv ll Ar 2 (i)- Subsection 15.21 explains how to introduce topology 
on this space. Let B(IL T ) be the corresponding Borel cr-algebra. Denote Pjy N ^ the probability 

measure on (n T ,£(n T )), induced by the process (Mzv lt jv 2 (i), < t < T). 

Our assumption for the theorem is that for any test function cp(x) the sequence {Mn 1) n 2 (0), 4>) 
weakly converges as Aq, N2 — > 00. 

We want to prove that as N± , N2 — > 00 the sequence of probability distributions Pjy Ni has 
a weak limit, and this limit is a one-point measure, that is the only trajectory (mi(t), 7712(f)), 
< t < T, which is the classical solution of the system ([TJ)— (J2J) - We split a proof of this result 
into the next two propositions. 

Proposition 6 The family of probability distributions |-P/Vi N2} on (^^' ^0^ T )) * s tight. 

Proposition 7 Limit points of the family of distributions Pjy N are concentrated on the weak 
solutions of the system (QP-fiP- 
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4.1.1 Tightness 

Before proving Proposition [6] we start with some preliminary lemmas. We want to prove that 
the family of distributions N2 of the random process (MN lf N 2 (i), < t < T), with values in 
the space of generalized functions, is tight. By the theorem 4.1 of [9j (see also Subsection 15. 2p . 
it is sufficient to prove that for any test function tjj = (ipi(x) , ^(x)) the family of random 
processes ((Afjv^JVaC*)* ifi) , < t < T), with values in R, , is tight. This will be done in the 
Proposition [10] below. 

Fix some test function ip = (ipi(x), ip2 (x)) and consider the random process 

^,JVi,JV a (a; (1) (t);x (2) (t)) = (M NltN2 (t),i(>) = 

= wX^ 1){t)) + wX Mx ? {t)) 
1=1 j=i 

This is a function of the Markov process £jvi jsr 2 (t), thus (see [6] Lemma 5.1, p. 330], for 
example) the following two processes are martingales: 



W A N liN2 {t) = F^, Ni ,n 2 (/'(^^(ijJ-F^^ (s (1) (0);x( 2 >(0) 

-J L NltN2 F^, tNuN2 (x (1) (s); x^(sfj ds (6) 

V^,^ N2 {t) = (W^ NuN2 (t)) 2 - J L NltN2 F^ NuN2 (x (1) (s);x (2) (s)) ds + 

+2^ F^ NuN2 ^'(s);^ 2 ^)) L Nl>N2 F^ NltN2 (a?«(*);aj( 2 >(*)) ds. 

For shortness we will write F(x^ ; x^) instead of -F^a^a^x^; x^). 
Lemma 8 The following estimates hold: 

i) \Lnx,n 2 F (x^; x^*) I < Ci("0, U2, Q!i2j «2l) uniformly in Ni, N% and (x^ 1 '; x^); 

ii) uniformly in x^\ x^ 

i^ w ;i w )-f(^ii w )i^ w ;i w 



Cl 2 («12, 4>l) C 2 l{a2l,1p2) 

(7) 



Proof of the lemma. Note that 



f((^)^-f(^) = i-(*(»f) -*(»!•>)), 

F((x«;x«) j )-F(x«;x«) - -*(*))• 



Thus 



ATi AT 2 

l Ki , w ,f _ * E « (,<■>) + « E «(»f) + 

1=1 j=l 
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iVi N 2 



i=l j=l 
AT 2 ATi 



+ 



+ 



EE^M^)-^))- 



i=i i=i 



(8) 



Then 



L^^Ffx^x^ 



< 



i Mc 



+ 



2||C« 



+ 



+2a i2 HV'ilIc + 2a 2i ll^llc 

and the assertion i) of the lemma is proved. To prove assertion ii) it is convenient to represent 
L N lt N 2 = L N!,N 2 +L Ni,N 2 as tlie sum of "differential" L^ ^and "jump" L 1 ^ ^ parts. 
It is easy to see that 

L° NltNa F°(zW;xW) - 2F( X U-, X M)L% uN2 F( X U-, X W) = 0. 

Let us prove that uniformly in x^ , x^ 2 ** 



tN2 F 2 (x« ; ) _ F(x« ■ x ( 2 ) ^ ^ F(x« ; ) 

In fact 



c 4«i2 H^illc + 4«2i \\ih\\c ^ 



No 



F 2 ^ X W- X (^^-F 2 (x«;x( 2 )) = ^(xtV^jJ-F^;^) 

x ( 2 F (x< V (2) ) + ^ (>i (xf)-Vi (x«)) 



Fl (x«;x (2) 



2F (x«;x( 2 ) 



i-tj 

2 



f(x^:x^ 



and similarly for expressions with (x^; x^). .. Thus 



Ni N 2 



"1 J»2 r / \ 

L^ N2 F\x^-x^) = 2F(x«;x( 2 ))^^^. F ( (*«; x^ .) - F (x«; x( 2 )) 



Nx N 2 

EE 

i=l 3=1 



x - x - Ol2 



i=l j=l 
1 



„(2) 



— ( i>i ( ) - ipi ( z 



(i) 



1 2 



+2f(x( 1 );x (2) )££^- F((x«;x( 2 )) 



j=l i=l 



N 2 N 



+EE^f- £W) -*(*?)) 



n^3 
1 2 



f(x«;x (2) 



j=i i=i 
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Nx f-f f-< A^ 2 iVi 



Af 2 Ni 



iV 2 f-f f-f iViAT 2 

the estimate ([9]) follows from this. Lemma is proved. 
Corollary 9 

sup E {W^ NuN2 (t)f -)■ 0, JVi, A 2 -»• oo. 

Proof. As P^i jVi JV 2 is a martingale with mean zero, it is sufficient to prove that the expectation 
of 



~L NuN2 F 2 (x«( S );x( 2 )( S )) - 2F (x^\s); x( 2 )( S )) La^F x( 2 )( S )) 

tends to zero. This follows from the estimate ([7]) of the lemma. 
Proposition 10 The sequence of distributions of real valued random processes 

Fi>,m,N 2 (x^(t);x^(t)) ,te[0,T], 



ds 



is tight. 

Proof of Proposition 1101 Remind that the following representation holds 

F^N U N 2 (x^(t);x^(t) 



F^, Nl ,N 2 (x (1) (0);x^(0)J + W^ NuN2 (t) + 
+ / L Nl>N2 F^ >Nl)N2 (x (1) (s);x (2) (s)) ds 



Note that our initial assumption is that the sequence -F^a^a^ (x W(0);x( 2 )(0)) weakly con- 
verges as Ni, N 2 — > 00. 

Prove now that the sequence 



n N ^(t) = j\ NuN2 F (x«( S );x( 2 )( S )) ds, t G [0, T] 



Ni,N 2 



is tight. We use subsection 6.1 of the Appendix. By assertion i) of the lemma 



L NuNa F(x^\s))x^\s)) ds 



< Cx(ip,v 1 ,v 2 ,ai2,a 21 ) ■ T , 



thus, the condition 1) of the Appendix holds. The condition 2) also holds, as one can prove 
that 

w'(rj Nl ' N2 ,-f) < 27 • C 1 (ip,vi,v 2 ,ai2,a 2 i) . 
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Prove that the sequence {W^^ 1 ,N 2 (t)y t G [Oj^]}^ j\r 2 * s tight. Using Kolmogorov's in- 
equality for submartingales with right continuous trajectories (see [3]), we have the following 
estimate, uniform in N\,N 2 , 

Then from the corollary [9] the condition 1) of Appendix holds. Thus 

P(\W (>r + ft\ W < E (W^ NuN2 {t + 9)- W^ NltN2 {r)) 2 



r+6»- 



E J r r Vf tNltN2 (s)ds 



e 2 



e 2 

Using this estimate one can check the sufficient condition of Aldous. Then Proposition 1101 is 
proved. 

This concludes also the proof of Proposition 
4.1.2 Weak solutions 

Definition 1 1 We say that the pair of functions M (t) = (mi (t , x) , m,2 {t,x)) is a weak solution 
of the system (EP - (EP? if for any pair 4>i(x), 4>2(x) £ S(R) the following identities hold 

(M{t),4>) = (M(O),0) + 

+ / (M(s), (ui0' 1 -Q!1201 + Q:2102,U202- "1201 + <*21 fa)) ds, 

Jo 

where 4>(x) = {<f>i{x) , ^{x)) , and the action of G(x) = (gi(x), g2{x)) on the test function 4>(x) 
can be written as 

(G,4>) = / gi(x)4>i(x) dx + / g2{x)4>2(x) dx 



Note that from the representation ([6]) and the identity j8|) it follows that 

{M NuN2 (t),(t>) = W t p jNl>N2 (t) + (M NuN2 (0),4>) + 

+ / ( M N 1 ,N 2 (s),{vi(/)[-ai2(pl + a2l4 , 2,V2^2-0!l24>l+a2l4>2))ds, 

Jo 

Let h = hit) G n T = D([0,T],S'(R) x S'(R)). For fixed <f> define the functional 

(h(t),(f)) - (h(0),(f)) - / (h(s), (viftx - a 12 (pi + a 2 i<p 2 , v 2 ^ 2 - 012^1 + a 2 i02))ds 



J<t>,r(h) = sup 

t<T 

In particular, 

SUplW^iVi.iVaWI = J((>,t(M Ni)N2 ). 
t<T 

The rest of the proof is standard (see [6]) and consists of three steps. 
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Step 1. From the definition of topology on II T it follows that J^tO) '• n T — > R + is a 
continuous functional. 

Step 2. Note that 

Ve > P{ J^HMnuN,) >e} = P£ uN2 {h : J^ T (h) > e} {Nt,N 2 -> oo) 
by Kolmogorov inequality and Corollary [9j 

Step 3. As J^t(-) is continuous, then the set {/i : JsTih) > 0} is open in It follows now 
that for any limiting point P^ of the family |-P/Vi iV 2 } we nave 

p L l h ■ J<t>A h ) > £ } < HmsupP^ N {h : J^ T ih) > e} . 

Ni,N 2 

That is, for any e > we have P^ {h : J^^ih) > e} = 0. In other words, all limiting points 
P^ of the family < P£ N > have support on the set {h : J^^rih) = 0}, which consists of 

weak solutions of {JTJ)— d2J) - 

This completes proof of Proposition 

The problem of uniqueness of the weak solution of JI])-([2]) is quite simple because the 
system (HJ)- ([2]) is linear. In the Subsection 14.21 we shall see that this system of first order 
differential equations has a unique classical solution which can be obtained in explicit way. 

4.2 Time asymptotics for the continuous model 

We prove here Theorem [3j 

Define the means (mass centrum) Oj(i) = / xm,i(t,x) dx and variance (momentum of 

From ([I])-® we get the following equations for the means 

a'i = v\ + a±2 (a 2 — a±) 
a'2 = V2 + «2i (oi — 02) 

It follows that equation for (12 (t) — a\ (t) is closed and has the following solution 

a 2 (t) - 01 (t) = V2 — ( 1 - e-^+^A + (a 2 (0) - 01 (0)) e -^ +a ^ t . 

Ct\2 + «21 \ ' 

Thus 

V% — V\ 

a 2 (t) — a\(t) -> [t ->• +00) 

OL\2 + «21 

and similarly 

— Oi(t) -> ; (t -)■ +00) 

as «12 + «21 



The equations for variances are 
d\ = a] 

d 2 = a 2 \ (di - d 2 ) + a 2 \ (ai(t) - a 2 (t)) 



d\ = a\2 (d 2 - d\) + «i2 (a 2 (t) - ai(t)) 2 
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Or, equivalently 



^ («2idi + 012^2) = 2ai 2 a 2 i (a 2 (i) - ai(t)) 2 

d 2 
— {d 2 -di) = - {ai2 + a 2 \) {d 2 - di) + (a 2 \ - ai 2 ) {a 2 {t) - ai(t)) 



From this we get 

d 2 (t) — d\(t) — > const 

and 



, 2 

v 2 — v\ \ a 2 \ — ct\ 2 



d\ 2 + «21 / «12 + «21 



— (a 2 idi + a 12 d 2 ) ->• 2a u a 2 i ■ 

dt V a i2 + a 2i/ 

Thus the growth of variances is asymptotically linear. Moreover, both are asymptotically 
equal. 

Now we come to the solution of the equations. Define the Fourier transforms 



i(x,t) = J exp(ixp)gi(p,t)dp 

We get 



— + Viipgi = a 12 {g 2 - gi) 



dg_ 
dt 

dg 2 . . , s 

— + v 2 ipg 2 = a 2 i{g! - g 2J 



with initial conditions m^O, x) = m i (x), i = 1,2. We write this system in the vector form 

dg , 
dt= A9 

where 

/ —ivip — a\ 2 a\ 2 



A 

For eigenvalues we have 



«2i —iv 2 p — a 2 ± 



a . la 2 



where 

a = i(v\ + v 2 )p + 012 + «2i, b = -v\v 2 p 2 + ip(y\a 2 \ + ^2012) 
One can write the solution as 

g = C + 4> + exp(tA + ) + C_0_ exp(iA_) 

where (f>± are eigenfunctions. Note that for small p there are two roots. One has ReA_ < 0, 
thus strongly decreasing term. Another is 

A+ = c lP + c 2P 2 + 0(p 3 ),c 2 / (10) 
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for small p. 

Let £t be a random variable with density m(x,t), g{k) - its characteristic function. We are 
interested in ^j(£t — a), a = E£t, its characteristic function is 

exp( " m ^ )5( 7! ) 

Using ([10p we get the result. 

Remark 12 One can see that there is no solution of the type 

mi(t,x) = fi(x - vt) 

as then fi would be exponents. 

Remark 13 For the singular initial conditions, that is when x, (0) = for k = 1, . . . , iV$; 
i = 1,2, one can get the same asymptotic results. 

5 Appendix 

5.1 Probability measures on the Skorohod space: tightness 

Let G [0, r])} ngN be a sequence of real random processes which trajectories are right- 

continuous and admit left-hand limits for every < t < T . We will consider £ n as random 
elements with values in the Skorohod space D T (R) := D ([0, T], R 1 ) with the standard topol- 
ogy. Denote Fj, the distribution of £ n , defined on the measurable space (Z?r(R), B (-Dr(R))). 
The following result can be found in [lj. 

Theorem 14 The sequence of probability measures {Pp} ng ]\j is tight iff the following two 
conditions hold: 

1 ) for any e > there is C{e) > such that 

supPp ( sup > C{e) J < e; 

n \0<t<T / 



2) for any e > 



lim lim sup P£ (£. : u/(£;7) > e) = 0. 

7^0 „ 



where for any function f : [0, T] — )■ R anc? any 7 > we define 

*"'(/; 7)= in | max sup |/(t)-/(s)|, 

moreover the inf is ewer a// partitions of the interval [0, T] suc/i i/tai 

= to < *i < • • • < U = T, U — ti-\ > 7, % = 1, . . . , r. 
The following theorem is known as the sufficient condition of Aldous [6j. 
Theorem 15 Condition 2) of the previous theorem follows from the following condition 

Ve > lim lim sup sup (|£ r+ - £ T | > e) = , 

n Ten T ,e<-y 

where TZt is the set of Markov moments (stopping times) not exceeding T . 



14 



5.2 Strong topology on the Skorohod space. Mitoma theorem 

Remind that Schwartz space S(R) is a Frechet space (complete locally convex space, the 
topology of which is generated by countable family of seminorms, that implies metrizability, 
see |llj). In the dual space 5"(R) of tempered distributions there are at least two ways to 
define topology (both not metrizable): 

1) weak topology on S"(R), where all functionals 



We shall consider <S"(R) as equipped with the strong topology Details can be found in 

The problem of introducing of the Skorohod topology on the space Dt(S') := D([0, T], S'(R)) 
was studied in [9] and |12| . The topology on this space is defined as follows. Let {pa} be a 
family of seminorms, which generates strong topology in S"(R). For each seminorm pa define 
a pseudometrics 



where the inf is over the set A = {A = A(i), t £ [0, T]} of all strictly increasing maps of the 
interval [0, T] into itself. Equipped with the topology of the projective limit for the family 
{oIa} the set Dt(S') becomes a completely regular topological space. 

Let B{Dt{S')) be the corresponding Borel a- algebra. Let {P n } be a sequence of probability 
measures on (Dt(S'),B(Dt(S'))). For each <j) £ S(R) consider a map : y £ Dt(S') —> 
y.(</>) £ Dt(R)- The following result belongs to I. Mitoma [9]. 



Theorem 16 Suppose that for any <j> £ S(R) the sequence ^Pnl^ j is tight in in Dt(R)- 
Then the sequence {P n } itself is tight in Dt(S'). 
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